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VIBRATION  OF  A  CANTILEVER  BEAM  THAT 
SLIDES  AXIALLY  IN  A  RIGID  FRICTIONLESS  HOLE 

INTRODUCTION 

Consider  a  cantilever  beam  that  fits  snugly  into  a  frictionless  hole, 
Fig.  1.   Axes  x,y  are  fixed.   A  force  F(t)  moves  the  beam  axial ly.   Hence, 
the  length  L  of  the  beam  outside  of  the  hole  is  a  function  of  time  t.   The 
lateral  deflection  of  the  beam  is  y(x,t),  where  for  x  <  0,  y  =  0.   The 
axial  velocity  of  the  beam  is  L  =  dL/dt. 

By  Hamilton's  principle, 


/ 


(5T  +  5W)  dt  =  0  (1) 


where  T  denotes  the  kinetic  energy  of  the  system  and  W  denotes  the  work 
of  external  and  internal  forces.   The  kinetic  energy  is 


T  =  j   pLQL2  +  i  p  /   yt2dx  (2) 


where  the  subscript  t  denotes  partial  differentiation,  p  =  the  mass  of  the 
beam  per  unit  length  =  constant,  and  L  =  the  length  of  the  beam  when  t  =  0 
The  strain  energy  of  the  beam  (the  negative  of  the  work  of  internal  forces 
for  a  conservative  system)  is 


/ 


U  =  I  El  /   y2   dx  (3) 

2    /    xx 


where  E  =  the  modulus  of  elasticity  of  the  beam  and  1  =  the  second  moment 
of  the  beam  cross-sectional  area.   Hence,  the  virtual  work  5W  is 

6W  =  F6L  -  6U  (k) 

where  6U  =  the  first  variation  of  the  strain  energy  (Eq.  3)-   Because  of 
the  inertia,  the  axial  movement  of  the  beam  tends  to  bend  the  beam  by  beam- 
column  action.   This  effect  is  neglected  here,  but  it  could  be  included. 
Equations  (1),  (2),  (3),  and  {k)    yield 


1 

/[pL  L6L  + 
° 


i  pyt2(L,t)  5L  +  F6L  -  1  Ely2x(L,t)  5l_ 


+  p  /   y  6y  dx  -  El  /   y   6y   dx]  dt  =  0 

/       t    t  /      XX   XX 


.5) 


If  L(t)  is  considered  as  given,  the  terms  with  the  factors  6l_  and  cL  merely 
determine  F  (_t )  to  give  the  prescribed  motion  !_(t).   The  Euler  equation  for 
the  lateral  motion  of  the  beam  is 


Ely     +  Py,v  =  0  (6) 

xxxx      tt 


The  natural  boundary  conditions  are 


»xx(L't)  =  W(L't)=0  (7) 


The  forced  boundary  conditions  are 

y(o,t)  =  y  (o,t)  =  0  (8) 


x 


Adopting  Hamilton's  viewpoint,  we  may  suppose  that  y(x,o)  and  y(x,t,) 
are   prescribed  (Fig.  2).   The  problem  is  to  determine  y(x,t)  in  region  R 
(Fig.  2)  subject  to  the  differential  equation  (Eq.  6)  and  the  indicated 
boundary  conditions  (Fig.  2). 

The  curved  boundary  at  the  right  side  of  region  R  is  troublesome.   It 
is  possible  to  transform  coordinates  so  that  this  boundary  becomes  straight 

Let 

x  =  x(£,n) ,  t  =  t(c,n)  (9) 

Regard  y  as  a  function  of  (£,r|).   Then, 

y  =  y(£,n),  5  =  ?Cx,t),  r,  =  n(x,t)  (10) 

and 


Wt  l3rnVt      Wc  dxjt 

dt  X  dt,    T)    dt  X     dH  E,  dt  X 

\  2jx  ^Vat'x   W^Wx 
d  t 

Here,  we  regard  y  as  being  a  function  of  (£»n,)- 
Now  let 


5  -  i  ,   t  =  n  0  2) 


When  we  write-^X  ,  -s~  ,  etc.,  we  imply  y  =  y(£,f|)  .   Hence,  -«—  =  (--£-)  (t-)  , 
and  s  imi 1 ar 1 y , 


3  y  _  1    3  y 

7~ IT  "  T   T 

3x    L    3? 


(13) 


Now, 


(h.)    =h.  (ii)    +  ly  =  y 

v3t;x       35    V3rx       3n       y 

and    since 

/3Cn               *L            5L        /3ri\ 
^3t^x   "         ,  2              L     '    V3t;x 

=  1, 


(H) 


we  find 


Y    d?;n  l   W? 


and  hence, 


35  n 


2  -2 

3^2;  l  v35;n  l   3C3n 

2    •  ••  -2 

3  y  ik  / £y\  5k  .  (§x.\    ik_ 

353n  L  v35;n  L   ^35;n  l 


(15) 


a2 
3n 

Thus,    by    the    third   of    Eqs.     (11),    Eqs.     (12),    Eqs.     (U),    and    Eqs .     (15),    we 
find 


(Lx.)    =  r-(i-X)  ik  .  (li)i  +  i_y_i  r.  ik) 

\    2;x        L    \_2;     L  ^35JL        3?3nJ  l       L  ' 
d  t                          at, 

L    k3£3r/  l  ^3£;  L  l3C  .  2        ,    2J 


L  3n 


and    since   n   =    t ,    we   obtain 


2       2   2*2         '2      2    • 

V  *  '  k2    l2  a?    l2        ^3t  L 


ay  c'l  ,  a2y 

*  L    3t2 


(16] 


in  which  y  =  y(^,t)  on  the  right  hand  side  of  Eq.  (16) 
Consequently,  Eqs.  (6),  (13)  and  (16)  yield 


k  2        2-9  -2 


^    V    l2      2*     l2 
-  9  l!x_  ik     9y  c'i.  ,  a2y-i 

3C8t     L     '    3U         .    2J 

0  t 


(17: 


where  y  is  now  regarded  as  a  function  of  c,   and  t.   Equations  (7)  and  (8) 
become 


yrr(l,t)  =  yrrrCt)  =  0 


y   (0,t)  =  y   (0,t)  =  0 


(13) 


Equation  (17)  is  much  more  complicated  than  Eq .  (6).   However,  the 
region  R  is  simpler  (Fig.  3)-   The  problem  remains  a  linear  boundary-value 
problem  if  in  Eq.  (17)  we  regard  L(t)  as  a  given  function  of  time  t.   If 
we  consider  L(t)  as  unknown,  we  must  obtain  another  equation  which  defines 
L(t).   This  equation  is  obtained  from  Eq .  (5)  from  the  terms  with  the 
factors  5l_  and  6lL.   Hence,  we  find 


L  +  2Fr[EI4a,t)    -py/(L.t)]-^F(t)  (19) 

O  O 


In    terms   of    (5»t),    Eq.     (19)    becomes 

L  +II-C^.y|5(i>t)-  p[-  ~yc(1,t)+  yt(l,t)]2} 
F(t) 


(20) 


pLo 


Equation  (20)  is  subject  to  the  initial  conditions 

L(0)  -  Lq 


(21) 


L(0)  =  0 

where  L   is  the  initial  length  of  the  beam  at  time  t  =  0,  (see  Fig.  1). 
Thus,  if  y(x,t)  and  L(t)  are  considered  as  the  unknowns,  Eqs.  (17) 
and  (.20)  subject  to  the  boundary  conditions  of  Eqs.  (18)  and  (21),  respec- 
tively, represent  a  complicated  nonlinear  boundary  value  problem. 

ANALYSIS  OF  EQUATIONS 

(a)  l_(.t)  Prescribed.   If  L(t)  is  given,  then  the  problem  reduces  to 
seeking  a  solution  of  Eq .  (17)  subject  to  the  boundary  conditions,  Eqs.  (18) 
The  problem  is  still  of  considerable  difficulty.   However,  the  problem  is 
a  linear  boundary-value  problem.   The  nature  of  the  solutions  of  Eq .  (17) 
is  unknown.   Indeed,  one  must  consider  whether  or  not  a  unique  solution  of 
the  boundary-value  problem  exists.   Also,  is  there  a  unique  solution  of 
the  corresponding  initial  value  problem  in  which  the  condition  y  =  F. (£) 
on  t  =  t.  is  discarded  and  the  functions  y(£,0)  and  y  (^,0)  are    prescribed? 
One  way  to  approach  such  questions  is  through  the  calculus  of  finite 


differences,  in  which  the  differential  equation  is  replaced  by  a  linear 
partial  difference  equation.   The  existence  of  a  solution  of  the  difference 
equation  is  a  question  of  linear  algebra.   The  most  practical  way  of 
developing  numerical  solutions  of  Eqs.  (17)  and  (18)  is  an  open  question. 
One  possible  approach  is  through  the  calculus  of  finite  differences,  or  a 
combination  of  the  calculus  of  finite  differences  with  the  Runge-Kutta 
method.   Another  approach  is  through  finite-element  methods. 

The  problem  is  simplified  greatly  if  the  axial  motion  is  so  slow 
that  time  derivatives  of  l_(t)  are  negligible.   For  example,  if  the  beam  is 
vibrating  in  a  natural  mode  and  it  is  then  drawn  slowly  into  the  hole, 
does  it  continue  to  vibrate  in  a  natural  mode?   How  do  the  frequency  and 
amplitude  vary  with  L  under  these  conditions?   For  a  rigid,  frictionless 
hole,  the  total  mechanical  energy  is  constant,  since  no  energy  is  imparted 
at  the  root.   This  criterion  serves  to  determine  how  the  amplitude  varies. 
However,  we  should  not  have  to  apply  this  criterion,  since  constancy  of 
mechanical  energy  should  follow  automatically  from  the  solution  of  Eq.  (17) • 

If  L  and  L  are    negligible,  Eq.  (J  71  reduces  to 

II  +  ^il=0  (22) 

3^   El   3t2 

Equation  (22)  is  much  simpler  than  Eq.  (17)-   However,  the  fact  that  L  is 
a  prescribed  function  of  t  still  complicates  it.   There  is  a  question  as 
to  whether  the  character  of  the  solution  is  altered  by  the  reduction  of 
Eq.  (17)  to  Eq.  (22).   Again,  does  the  initial  value  problem  posed  by 
Eqs.  (18)  and  (22)  have  a  unique  solution?   If  so,  the  constancy  of  mechanica 
energy,  as  a  separate  principle,  is  not  needed  to  determine  the  dependency 
of  amplitude  upon  L.   In  fact,  constancy  of  mechanical  energy  should  be  a 
deducible  consequence  of  Eqs.  (18)  and  (22). 


Unlike  Eq .  (17),  Eq.  (22)  allows  a  theory  of  natural  modes.   For 

example,  consider  solutions  of  Eq.  (22)  of  the  form  y  =  f(£)g(t).   Then, 
Eq.  (22)  yields 

f""  ,  pl!  g^_  Q 

f  El   g 


Hence, 


L 

fun    u      ni  n   nii      h 

f-=s-trf-=-s  (23) 


in  which  6  is  a  constant.   Then, 


g..  +  e_Eig  =  0  (24) 

PL 


Unfortunately,  Eq.  (24)  does  not  yield  simple  harmonic  motion,  because  L 
is  a  function  of  t.   However,  the  equation 

f""  -  B  f  =  0  (25) 

can  be  integrated.   With  the  root  conditions  (Eqs.  18)  f(o)  =  f'(o)  =  0, 
Eq.  (25)  yields 

f(g)  =  A(sinh6?  -  sin6?)  +  B(cosh6C  -  cosBC)  (26) 

Now,  the  free-end  conditions  (Eqs.  18)  f"(l)  =  f"r(l)  =  0  yield 


A(sinhB  +  sinBL  +  B(cosh6  +  cosBl  =  0 
A(coshg  +  cos3L  +  B(_sinh6  -  sing)  =  0 

and  the  vanishing  of  the  determinant  of  Eqs.  (27)  yields 


(27) 


coshB  cos3  =  -1  (28) 


These  results  are  exactly  like  that  for  natural  modes  of  a  cantilever  with 

fixed  length.   The  roots  of  Eq.  (28)  are   constants  B,  ,  B~  ,  B,  ,  •••  .   They 
have  been  tabulated  by  Young  and  Felgar  (1). 

As  is  customary,  let  B  =  1  and  A  =  -a  .   Then,  by  Eq.  (27), 

coshB  +  cosB 

a  =  - -  (29) 

n   sinhB  +  sinB 
n       n 

Consequently,  a.,  a_ ,  •••  are   constants,  independent  of  L.   The  values  of 
a  have  been  tabulated  by  Young  and  Felgar.   Equation  (26)  gives  the  natura 
modes , 


f  (?)  =  coshB  5  -  cosB  C  -  a  (sinhB  5  -  sinB  O  (30) 

n  n        n     n      n        n 

The  functions  f  {£,)    are   orthogonal  in  the  interval  (0,1);  in  fact,  it  is 
eas  i 1 y  shown  that 


(3D 


The  fact  that  the  functions  f  (E)  satisfy  the   end  conditions  (Eq.  13), 

n 

makes  then  ideal  approximating  functions  for  use  in  numerical  attacks  on 
Eqs.  (17)  and  (22).   Presumably,  a  solution  of  Eqs.  (17)  or  (22)  can  be 
approximated  by 


yC5,t)  =  ]T   qn(t)fnU)  (32) 

n=l 

This  approach  reduces  the  problem  to  one  of  a  finite  number  of  degrees  of 
freedom.   The  generalized  coordinates  are   q  (t).   The  Lagrange  theory  of 
linear  vfbrations  is  then  applicable.   Accordingly,  let 

6  2   — 

u>  --4-/U.  (33) 

n    L2     p 

Then,  Eqs.  (24)  becomes 

g  "  +  a)  2g  =0  (3^) 

3n     n  3n 

Since  L  depends  on  t,  note  that  oo  depends  on  t.   Consequently,  co   is  not 
exactly  the  frequency  of  a  harmonic  motion.   The  general  solution  of  Eq. 
(3M  is 


g  =  g  (t,  A  ,  B  )  (35) 

3n   3n     n   n 


in  which  A  ,  B  are   constants  of  integration.   An  infinite  series  solution 
n   n 


of  Eq.  (22)  is  then 


Zf  (C)g(t,  A,  B)  (36) 

n    n     n   n 


n=l 


Equation  (36)  automatically  satisfies  boundary  conditions  (Eqs.  18)  at  the 

free  end  and  the  root.   One  might  expect  that  the  constants  (A  ,  B  )  can  be 

n    n 

chosen  to  satisfy  the  conditions  at  t  =  0  and  t  =  t .;  i.e., 


10 


V*  f  (C)gn(0,  A   B  )  =  Fn(0 

^_^   n    n     n    n     0 


n=1 


(37) 


z 

n=l 


'n^W  VBn] 


=  F,(?) 


Alternatively,  the  constants  A  ,  B  might  be  chosen  so  that  the  beam  has 

n   n   3 

a  given  initial  deformation  (the  first  of  Eqs.  37)  and  given  initial 
velocities;  i.e.,  y  (5,0)  might  be  specified  instead  of  y(£,t,). 

If  L(t)  is  a  sufficiently  slowly  varying  function  of  t,  the  approxi- 
mate  solution  of  Eq.  (34)  is 


g  =A   sincot  +  B   coso:  t 
n    n     n     n     n 


(38) 


Then,  Eq.  (.361  becomes 


y(5,t)  =  7    (A   sinw  t  +  B   cosoj  t)  f  (5! 
'        /    j  n     n     n     n    n 


n=l 


The  initial  value  problem  then  reduces  to 


(39) 


n-l 


Bnfn(-} 

n  n 


-  Fo(5) 


Co) 


n-l 


VA^)  "  +o(5) 


in  which  it  (£)  =  y  (£,0),  and  oj    is  the  initial  value  of  uj  .   Under  fairly 
o       t  n  n 

broad  conditions,  the  constants  A  ,  B   can  be  chosen  to  satisfy  Eq .  (kO) 

n    n 

when  F  (£;)  and  (j)  (£)  are  prescribed  functions.   Then  Eq.  (39)  represents  the 
solution  of  the  initial  value  problem. 


The  kinetic  energy  and  the  strain  energy  are,  respectively, 


00 

-f      ^=!pLZ  ^n 

i  n=1 


T   =  7T  pL  /        Y^dE,    =  -^  pL   y         (jo   ^(A      coseo    t    -   B      since   t) 
2         /  t  2  /   j        n        n  n  n  n 


1 


El/  2  El  V    >  h 

U   =  -=-=-  /        yrrd5   =  — =-      >        8      (A      sinco   t   +   B      cosco   t 
o,  3    /  KK  n,  3      fc    ^       n        n  n  n  n 


21/  /     ^     2L 
0 


Hence,  in  view  of  Eq.  (33),  the  total  mechanical  energy  is 


T  +  U  =  —-V  6  \a  2  +  B 
2,3  L   4        n    n     n 


2) 


n=l 
Thus,  the  total  mechanical  energy  is  independent  of  t,  provided  that  A 


n 


3/2 
and  B  vary  as  L    .   Apparently,  for  the  conditions  stated,  this  require- 
ment ensures  conservation  of  mechanical  energy. 

(b)  F(t)  Prescribed.   If  one  prescribes  F(t),  then  Eqs.  (17)  and  (20) 
subject  to  Eqs.  (18)  and  (21)  must  be  solved  simultaneously  in  order  to 
determine  l_(t).   The  method  of  solution  for  both  cases  l_(t)  prescribed 
and  F(t)  prescribed  is  outlined  in  the  following  section. 

Method  of  Solution 

Equations  (17)  and  (20),  together  with  the  initial  and  boundary  condi- 
tions given  by  Eqs.  (18)  and  (21)  form  a  system  of  two  nonlinear  partial 
differential  equations  in  the  two  unknowns  l_(t)  and  y(x,t).   An  approximate 
numerical  solution  of  these  equations  can  be  obtained  by  a  Galerkin  finite 
element  formulation. 


1  2 


Thus,  in  accordance  with  the  Galerkin  FEM,  we  form  the  approximate 
solutions  of  L  and  y  as  $  and  v  which  are, 

L(t)  =  <J>(t)  =  8(t) 

T  (M) 

y(x,t)  =  v(x,t)  =  N  (x)a(t) 

where  N  are   a    set  of  quadratic  shape  functions  with  local  support.   A 
quadratic  shape  function  is  associated  with  each  nodal  point  of  the  dis- 
cretized  finite  element  model.   Quadratic  interpolation  over  an  element 
requires  each  element  has  three  nodal  points,  say  one  at  each  end  of  the 
element,  and  a  nodal  point  in  the  center  of  each  element. 

Figure  4  shows  that  the  odd  numbered  shape  functions  N.,  N_ ,  ..., 
associated  with  nodal  points  at  the  ends  of  an  element  span  five  nodal 
points  (two  elements),  while  the  even  numbered  shape  functions  N  ,  N,,  ... 
span  the  three  nodal  points  of  one  element. 

In  accordance  with  the  Galerkin  FEM,  we  form  the  residual  function 
for  Eqs.  (171  and  (20)  as, 


•2        -2 
RjCx.t)  =  ^  v""  +  p[v"C2^2  +  2v'S  4_ 


-  2v'?  £  -  v'5  ■£  +  v]  (42) 


R20.t)  =  i   +  2HT  [T"  v"2(1't)  "  p{"1  V(ift)  +  v(i,t)}2] 

o  4> 


^-F(.)  <*3> 

o 


The  finite  element  equations  are   obtained  by  requiring  that  the  residual 
functions,  given  by  Eqs.  (42)  and  (43),  be  orthogonal  to  each  of  the  basis 
functions.   That  is, 
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/ 


NR,dx  =  0 


(Mt) 


/ 


HR„dx  =  0 


(45) 


where  H  is  the  function  given  by 


H  =  1 


0  <  x  <  1 


(46) 


Substitution  of  Eq.  (4l)  into  (42)  and  (43)  gives 


•2  -2 

R,  (x,t)  =  A-(EIN"Ta)"  +  p[K2  K   N"Ta  +   2£  K   N,Ta 


'1 


TV 


2  "  ~    *   2  2     ~ 
6  B 


25  -I  N|Ta  -  5  I  N,Ta  +  NTa] 
p  ~       p  ~      ~~ 


(47) 


R,(I,t)  -  5  +  j±- [%N"Ta)2 


5=1 


-  p(-|N'Ta 


+  NTd|    }2] 


t,=  i 


PL 


F(t) 


(w: 


Substitution  of  Eqs.  (47)  and  (48)  into  Eqs.  (44)  and  (45)  respectively, 
and  performing  integration  by  parts,  gives 


1 


1 


1 


BT  +  /   EIN"N"Ta5a-  ^y  /  p?2N ' N ,Ta^a  +   2  —  j     ?NN|TdCa 
J  ~     B  J  ~  ~    ~     e^  /   ~~ 


0 


1 


-  2  I  /   ?NN|Td5d  -  |  /   ?NN|TdCa  +  /   NNTd?a  =  0 


(49) 
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and 


a  1   El  r /  ^  8         k  ,2 

2pLo  B»  I   2  2n_1   X,  2  2n   £  2  2n_1 
e  e         e 

-  p{-|(i-a2n_i  -  iLa^  +f  a2n+1)  +  a2n+1}2]  -  ^  F(t) 


(50) 


where  BT  denotes  a  boundary  term. 
Lett  ing , 

1 
A  =  /   El  N"N"T  d£ 


B  = 


p£2N'N|T  dE, 


C  = 


C  N  N,T  d? 


(50 


D  =  /      N  N   d£ 


3  =  (72"  Gri- 
lle 


8         4       ,2 

a_  +  — *■  a_  ,  . ) 


,  2  2n   n    2  2n+ 

£        £e 


b  =  {f  (f  a2n-l  "  r  a2n  +  T   a2n+l}  +  d2n+l}2 
e         e        e 
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Equations  (k3)    and  (50  become 


B         B 


BT  +  Aa  -  £=-  Aa  +  2  ^  Ca  -  2  f-  Ca  *  Da  =  0 


B 

6"  ~r 


(52) 


o  B  o 


(53! 


Defining  the  following  terms, 


Y  =  -(BT  +  Aa) 


a)  = 5-B  Aa  +  2-^-Ca 

B     ~     B 

C   =  -2  —  C 


(54) 


5  =  -  B  c? 


4-  F(t)  -  ~-   (a  -  pb) 


pL    v  '    2pL 


jhere  B   is  the  value  of  B  at  the  previous  time,  Eqs.  (52)  and  (53)  become, 


Da  +  eg  +  C  a  +  oaB  =  y 

B  =  G 

or,  in  matrix  form 

1         2n+l  2n+2 


2n+l 
2n+2 


(55! 


U) 


o  :  o 


(56) 


Equations  (55)  (or  (56)),  subject  to  initial  conditions  may  be 
solved  by  numerical  integration. 
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F(t) 


Figure  1.   Cantilever  Beam  Moving  Axially  in  Cylindrical  Hole, 


U=^=o 


1 


Region  ft 

D>FFEP.£NT!AL  E^.  (6) 


^V 


X) 


X=L(t) 


*,»=!)««"  ° 


Figure  2.   Region  of  Integration  for  Equation  (6) 
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t 


i  3=f.(s)  '• 


Figure  3.   Region  of  Integration  for  Equation  (17) 
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Figure  k.      Discretized  FEM  Model  with  Quadratic  Basis  Functions 
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